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ABSTRACT 

It was suggested that some of the short-duration Gamma-Ray Bursts (GRB) are giant flares of Soft Gamma- 
ray Repeaters (SGR) in nearby galaxies. To test this hypothesis, I have constructed a sample of 47 short GRBs, 
detected by the Inter-Planetary Network (IPN), for which the position is constrained by at least one annulus on 
the celestial sphere. For each burst, I have checked whether its IPN 3-C7 error region coincides with the apparent 
disk of one of 316 bright, star-forming galaxies found within 20 Mpc. I find a single match of GRB 000420B 
with M74, which could, however, be due to a chance coincidence. I estimate the IPN efficiency as a function 
of fluence and derive the galaxy sample completeness. I find that assuming there is a cut-off in the observed 
energy distribution of SGR flares at < lO'*^ erg, the fraction of SGRs among short GRBs with fluence above 

10^^ erg cm^^ is < 16% (< 27%) at the 95% (99.73%) confidence level. I estimate the number of active 
SGRs in each one of the galaxies in the sample, and combine it with the distances to these galaxies, the IPN 
efficiency, and the SGR flare energy distribution (Cheng et al.), to derive the rate of giant flares with energy 
above 4 x 10^^ erg. I find that the rate of such giant flares is about (0.4 — 5) x 10^^ yr^' per SGR. This rate 
is marginally consistent with the observed Galactic rate, of a single giant flare with energy above 4 x lO'*^ erg 
in 30 years. Comparison of the Galactic rate with the inferred extragalactic rate implies a gradual cut-off (or 
steepening) of the flare energy distribution at < 3 x lO'*^ erg (95% confidence). Using the Galactic SGR flare 
rate, I set a lower limit of one percent on the fraction of SGR flares among short GRBs. 

Subject headings: gamma rays: bursts: individual: (GRB 000420B, GRB 000526B, GRB 051103) — stars: 
neutron — galaxies: individual (M74, NGC 7331) 
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1. INTRODUCTION 

The rate of Soft Gamma-ray Repeater (SGR) giant flares 
and their fraction among short-duration Gamma Ray Bursts 
(GRB; Kouveliotou et al. 1993) are important ingredients for 
the understanding of giant flares in the context of the magnetar 
model (Duncan & Thompson 1992; Paczynski 1992; for a 
recent review see Woods & Thompson 2006). 

The large energy release from the 2004 December 27 giant 
flare, combined with some similarities between the temporal 
and spectral properties of giant flares and short GRBs, has re- 
ignited the idea that some or all of the short GRBs are in fact 
SGR giant flares in nearby galaxies (e.g., Dar 2005; Hurley et 
al. 2005; Palmer et al. 2005; Nakar et al. 2006). Additional 
support to this hypothesis came from the observed Galactic 
rate of giant flares with energy above 4 x 10^^ erg, of one in 
30 years among four active Galactic SGRs. 

A limit on the fraction of SGRs among short GRBs was ob- 
tained by Nakar et al. (2006), who searched the error quadri- 
lateral of six, well-localized, Inter-Planetary Network (IPN) 
GRBs for relatively nearby galaxies. They put a upper limit 
of < 40% (95% confidence limit; CL) on the fraction of SGRs 
among bright short GRBs. 

Lazzati, Ghirlanda, & Ghisellini (2005), inspected the spec- 
tra of 76 Burst And Transient Source Experiment (BATSE) 
short bursts. They found three GRBs whose spectra are well 
described by a black body model, and 15 bursts for which a fit 
with a black body model is not excluded. Lazzati et al. (2005) 
then argued, based on the gamma-ray light curves of these 
short GRBs, that none of them were likely SGR bursts. They 
derived a limit on the fraction of SGR bursts among bright 
BATSE events of < 4% (95% CL). I note that the Lazzati et 

' Division of Physics, Mathematics and Astronomy, California Institute of 
Technology, Pasadena, California 91125, USA 



al. (2005) results are based on the assumption that all the gi- 
ant flares have similar light curves. However, this assumption 
is based on a small number of observed giant flares. If I relax 
the Lazzati et al. (2005) assumption about the temporal char- 
acteristics of giant flares, but keep their assumption regarding 
the spectral properties of SGRs (and there are up to 18 SGR 
flares in their sample), their result poses an upper limit on the 
fraction of SGRs among bright short GRBs of < 35%, at the 
95% CL. 

Recently, two possible detections of extragalactic SGRs 
were claimed: Crider (2006) detected a 13.8 s periodicity 
in the tail of the gamma-ray light curve of GRB 970110. 
He suggested that this burst may be an SGR flare in 
NGC 6946. Golenetskii et al. (2005) reported the detection 
of GRB 05 1 103, which includes the nearby galaxy M8 1 in its 
error quadrilateral (see Ofek et al. 2006 and Frederiks et al. 
2006 for follow-up observations and discussion). 

Another quantity related to the fraction of SGRs among 
short GRBs is the rate of SGR giant flares. This rate is the 
consequence of the total energy available for SGRs, which in 
the context of the magnetar model is provided by the magnetic 
field of neutron star 

Recently, Stella et al. (2005) argued that flie rate of SGR gi- 
ant flares, with energy above 5 x 10^^ erg, is about 10^^ yr^'. 
They obtained this result by using a uniform prior (see how- 
ever, Duncan 2001) on the observed rate. They further argued 
that about 70 giant flares, with energy above 5 x 10^^ erg, are 
expected during the lifetime of an SGR, and concluded that 
magnetic fields of B ^ lO'*' G are needed in order to explain 
the energy source of SGRs. 

Popov & Stern (2005) attempted to find extragalactic SGRs 
by searching the BATSE catalog for short GRBs that spatially 
coincide with the Virgo cluster. They used this search to put 
an upper limit on the rate of giant flares with energy above 
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Fig. 1. — Histogram of BATSE GRBs integrated four-channel fluence, 
for short-duration bursts triggered by BATSE (empty bars), and BATSE- 
triggered short-duration GRBs that were detected by at least two IPN de- 
tectors and thus constrained by at least one annulus (filled bars). 



~ 10'*'' erg, of 10^^ yr^'. However, this limit was obtained by 
assuming that Tso of all SGR giant flares is in the range of 0.05 
to 0.7 s, and by ignoring the BATSE positional uncertainty. A 
similar limit was obtained by Palmer et al. (2005). 

In this paper 1 match the IPN short GRBs with nearby 
galaxies (Sj3)- Contrary to previous efforts, I do not introduce 
any assumptions about the temporal or spectral properties of 
SGR giant flares. 1 present, in 2] the search for IPN short- 
duration GRBs that spatially coincide with nearby galaxies. 
In SjS] I use the results of this search to calculate the rate of 
SGR giant flares and confront it with the rate estimated based 
on the Galactic SGRs. 1 finally discuss the results in S|6] 

2. THE IPN SENSITIVITY 

The Inter-Planetary Network (IPN; e.g., Laros et al. 1997; 
Hurley et al. 1999) is a set of gamma-ray detectors on board 
several spacecraft in the solar system. The IPN uses GRB 
photon arrival times to triangulate their celestial positions. 
Each pair of spacecraft constrain the GRB position to an annu- 
lus on the celestial sphere. The area of IPN 3-a error regions 
is typically two orders of magnitude smaller than BATSE 
error circles. In fact, IPN positions were used to calibrate 
BATSE positional errors (Briggs et al. 1999). 

Here I have use the IPN catalog^ containing all IPN- 
triangulated GRBs observed from 1990 November 12, to 
2005 October 31^. Figure [T] shows a histogram of the 
integrated four-channel fluence for BATSE-triggered short- 
duration GRBs (Tgo < 2 s; empty bars), as well as all short- 
duration GRBs triggered both by BATSE and an additional 
IPN detector (i.e., that have at least one annulus constraint; 
filled bars). 

Figure |2] shows the approximate all-sky IPN complete- 
ness as a function of BATSE integrated four-channel fluence. 
Completeness here is calculated from the ratio of the two his- 
tograms in Fig. [T] multiplied by the exposure completeness 

^ http://www.ssl.berkeley.edu/ipn3/interpla.html; version 2005 De- 
cember 17 - For reference, this IPN catalog version is available from 
htpp://astro.caltech.edu/~eran/GRB/IPN/NearbyGaVCatIPN.txt.verl7122005 

^ This catalog does not contain GRB 051103 which may have originated 
in the nearby galaxy M81 (Ofek et al. 2006; Frederiks et al. 2006). 
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Fig. 2. — The all-sky IPN approximate completeness as calculated from the 
ratio of the two histograms in Fig.[T] multiplied by the exposure completeness 
(86%; see text). The upper x-axis shows the corresponding distance to a 
3.7 X lO''* erg event, which is the energy of the 2004 December 27 flare 
assuming a distance of 15 kpc to SGR 1806—20 (Hurley et al. 2005). The 
dashed line is an approximation to this completeness function (Eq.[T). 

of the KonusAVind and BATSE short-duration GRB sample 
that is presented in the next section. I estimated the exposure 
completeness to be about 86% (= [4 x 0.48 + 11 x l.Oj/15). 
This is based on the fact Konus/Wind started to work only in 
1994 (11 yrs out of 15 yrs), and the exposure completeness 
of BATSE, which worked 4 yrs without KonusAVind, is about 
48%^. The dashed line shows an approximation to this com- 
pleteness function given by: 

0, F<5xlO-^ 
0.3281ogioF + 2.065, 5 x 10-^<F<2.1 x lO""* 
0.86, F>2.1xl0-^ 

(1) 

where F is the fluence. BATSE sensitivity is almost complete 
down to a fluence of about 10^^ erg cm^^, in which the IPN 
efficiency drops to about ten percent. Therefore, the sensitiv- 
ity function (Eq. [T]i is an approximation to the IPN absolute 
mean efficiency. However, I will show in [J3] that the results 
are not very sensitive to the exact shape of this completeness 
function. 

3. SELECTION OF THE GRB AND GALAXY SAMPLES 

I selected from the IPN catalog all GRBs with a 3-(7 error 
region constrained by at least one annulus with semi-width 
smaller than one degree (1260). I further selected from this 
list all GRBs that were triggered by BATSE (752) and have 
Tgo < 2.0 s (29 events). The Tgo duration at which a GRB has 
an equal probability of being a short- or long-duration GRB 
is about 790 = 5 s (Donaghy et al. 2006). As SGR giant flares 
have durations shorter than 2 s, the upper limit on the fraction 
of SGR giant flares among short-duration GRBs derived in 
this paper is conservative. 

In addition, I searched for short-duration GRBs among the 
Konus/Wind satellite observations (Cline et al. 2003). A cat- 
alog of short GRBs detected by Konus/Wind was published 
by Mazets et al. (2004), however, some of the GRBs in this 

^ Calculated from the BATSE all-sky exposure map 
|http://cossc.gsfc.nasa.gov/docs/cgro/cossc/batse/4Bcatalog/4b-exposure.htrnll 
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Fig. 3. — The BATSE Tgo vs. the estimated duration (by eye from 
KonusA¥ind light curve), Tgy^, of the GRBs. 

catalog were classified as long GRBs by BATSE. Therefore, 
instead of using this catalog I used the following statistical 
approach: I inspected by eye the Konus/Wind light curves for 
which an IPN annulus constraint is available, and selected all 
the bursts that have a total duration of less than two seconds. 
Due to the differences between the BATSE and Konus/Wind 
detectors, this sample is contaminated by long GRBs (which 
appear like short GRBs for the Konus/Wind detector). Fig- 
ure [3] shows the BATSE Tgo vs. the duration estimated by 
eye (Tgyg) of GRBs that were triggered by both BATSE and 
KonusAVind. I found 83 Konus/Wind bursts with Tgye < 2 s 
and annulus semi-width of less than one deg. Of these 83 
events, 37 were detected by BATSE, of which 23 were classi- 
fied by BATSE as long-duration bursts (i.e., Tgo > 2 s). There- 
fore, the Konus/Wind sample has ^ 60% (23 /37) contamina- 
tion by long GRBs. Given this contamination factor, I esti- 
mate that the actual number of short bursts, among these 46 
remaining events, is 18^4. Finally, I compiled a list of 75 
events: 29 BATSE events with Tgo < 2 s and 46 Konus/Wind 
events with Teye < 2 s. This list contains about 47 (~ 18 -I- 29) 
short bursts. 

Next, I selected a sample of bright star-forming galaxies 
within 20 Mpc from a modified version of the Tully (1988) 
nearby galaxies catalog. The selected galaxies have absolute 
magnitude Mb < — 18, reside at distances less than 20 Mpc, 
and have morphological types other than E or SO. Introducing 
the last criterion, I assume that SGRs form predominantly in 
young stellar populations (see however, Levan et al. 2006). 
I excluded from this list the Large Magellanic Cloud (LMC). 
The nearby galaxy sample contains 316 galaxies, which are 
listed along with their basic properties in Table [T] Also listed 
in this table are the far Infra-Red (IR) fluxes and estimated 
supernova (SN) rates in these galaxies derived in [JS] 

3.1. The galaxy sample completeness 

There are three factors which affect the completeness of this 
galaxy sample: (1) galaxies behind the Galactic plane; (2) 
star formation in galaxies not included in the sample; and (3) 
the fact that the sample is volume limited. I estimated these 
incompleteness factors below. 

To estimate the completeness of the galaxy sample (Table[T]l 
due to Galactic obscuration, I calculated the number of galax- 
ies in this catalog as a function of Galactic latitude, b. I then 
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Fig, 4. — The survey volume completeness, Ey (Eq,|2j, as a function of 
distance r,„ov, for four cases: E„,ax = lO"" erg, a = 5/3 (solid line); E„,ax = 
5x10"^ erg, a = 5/3 (dashed line); E,„ax = 10"" erg, a = 9/5 (dashed-dotted 
line); and = lO*' erg, a = 5/3, and a modified completeness function 
(Eq,[T] C = C(0,8F); dotted line). 

extrapolated the number of galaxies from the unobscured re- 
gions (i.e., b ^ 30 deg) to the obscured regions. Based on 
this extrapolation, I estimated the completeness of the galaxy 
sample to be about eg = 90%. 

Next, I estimated, Esfr, the fraction of star formation within 
galaxies with Mb < — 18 relative to the total star formation 
within 20 Mpc. I calculated the far Infra-Red (IR) flux^ of 
the galaxies, defined by 2.58/60 + /10O' where /go and /loo are 
the 60 /xm and 100 /xm fluxes, respectively. This quantity is a 
good estimator of the star formation rate in galaxies (Helou et 
al. 1988). Given the galaxy distances, I summed the absolute 
far-IR flux of all the galaxies (excluding EO/SO) brighter than 
Mb < — 18 within ten Mpc and divided it by the absolute far- 
IR flux of all galaxies (excluding EO/SO) brighter than Mb < 
— 12. Note that I use a distance cut of ten Mpc in order to 
avoid incompleteness due to missing faint galaxies. I found 
that EsFR = 95%. 

The galaxy sample contains galaxies within 20 Mpc. How- 
ever, SGR giant flares can be observed to larger distances. 
To estimate the survey volume completeness, I divide the ex- 
pected number of SGR flares (assuming constant number den- 
sity of SGRs) that can be observed by the IPN up to a distance 
of r,„ax, by the total observable number of IPN SGR flares: 
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where r is the distance, C{F) is the completeness function 
(Eq. [TJ, and Emax is the break in the flare energy distribution. 
I weighted the number of flares of energy E, dN/dE, by the 
empirical E^" factor, with a = 5/3 (Cheng et al. 1996). 

FigurelHshows ey as a function of the distance r^ax^ for four 
different combinations of E^ax, OC and C{F) (see legend). It 
is clear that the volume completeness is not sensitive to small 
changes in a, nor to the exact shape of the completeness func- 
tion. However, the completeness, £v{rmax — 20 Mpc) is sen- 

^ Based on Infra-Red Astronomical SatelUte (IRAS) point source catalog, 
version 2,0, IPAC (1986) 
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TABLE 1 
Nearby bright galaxies 



Name 
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SN rate^ 
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0.00026 
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5 
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00 
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2.7 
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3 


7.853 


0.01685 
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01 


:36 


42 +15 


47 


00 


9.61 


-20.47 


7.0 


10.4 


5 


19.736 


0.01331 


N7331 


22 


:37 


06+34 


26 


00 


9.67 


-21.26 
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15.3 


4 


131.932 
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Note. — The first ten lines of TablefTI For convenience, I list also M74 (NGC 626) and NGC 7331 which are mentioned in this paper The entire table is available via the electronic 
version. 

^ Apparent blue magnitude, fij, corrected for Galactic extinction, and for the external galaxy extinction (see TuUy 1988).^ Absolute blue magnitude as calculated from the corrected 
distance and apparent extinction-conected magnitude.'' Distances are based on recession velocities, assuming Hq — 70 km s"^ Mpc"' The distances listed in the Tully (1988) catalog 
are based on galaxy recession velocities, assuming an Hubble constant Hq — 75 km s"' Mpc"'. I corrected all the distances in this catalog by a factor of (75/70; i.e., to conveit to 
Hq — 70 km s"' Mpc"' ), and recalculated the absolute magnitude accordingly."^ Galaxy type, see Tully (1988) for details.^ The total observed far IR flux, 2.58 /go + /loO' calculated by 
summing the fluxes of the IRAS sources found within the galaxy radius from the galaxy position.^ SN rate estimated by normalizing the absolute far IR flux to that of NGC 253, and 
assuming this galaxy has 0.2 SN yr"' (Pietsch et al. 2001). 




Fig. 5. — The survey volume completeness. By (Eq.[2), for r„„„ = 20 Mpc, 
as a function of the break in the flare energy distribution, £,„av (assuming 
a = 5/3). 

sitive to the break energy, E,„ax- Fig|5]shows the volume com- 
pleteness, Evirmax = 20 Mpc), as a function of E„,ax- A rea- 
sonable value for is 10"*^ erg, which roughly corresponds 
to the total magnetic energy available for a B = lO'^ G neutron 
star Assuming Ei„ax = lO'*^ erg, the volume-completeness 
of our survey is £v{rmax = 20 Mpc) = 0.75, while if E,„ax = 
5 X 10^^ erg, then the survey completeness drops to 30% (see 
however SjSj. 

4. THE SEARCH FOR SGR FLARES FROM NEARBY 
GALAXIES 

For each burst in the short GRB sample, I searched for over- 
lap between its IPN 3-(J error region and the apparent disk of 
one of the galaxies in Table [T] The IPN constraints include 
one or more annuli on the celestial sphere. In addition to 
the annuli, the IPN database occasionally contains constraints 
from individual instruments: ecliptic latitude constraints from 
Konus/Wind; error radii from various satellites; and/or planet 
blocking regions (i.e., in cases that Earth or Mars blocked part 
of the detector field of view). Some of these error regions 
(e.g., BATSE; WATCH Sazonov et al. 1998; COMPTEL Kip- 
pen et al. 1998; EGRET Gonzalez et al. 2004; PHEBUS 
Tkachenko et al. 2002; SIGMA Claret et al. 1994; Bep- 



poSAX Guidorzi et al. 2004; HETE Vanderspek et al. 2003) 
are given at the \-a confidence. For these error radii I added 
known systematic errors (i.e., 1.6 deg for BATSE; Briggs et 
al. 1999) and translated all the errors, which I assumed are 
circular^ and normally distributed, to 3-(7 errors. 

To avoid incompleteness due to the positional inaccuracies 
of the Tully (1988) catalog, the galaxy diameter used in the 
search was increased by one arcminute. 

I have found two matches between the short GRB can- 
didates and the galaxies listed in Table [T] The first, 
GRB 000420B, coincides with the position of the nearby Sc- 
type galaxy M74. Fig |6] shows a Galaxy Evolution Explorer 
(GALEX) near UV-band image, emphasizing the star-forming 
regions in M74, with the IPN annulus overlayed (solid lines). 
GRB 000420B was detected by Konus/Wind at UTC 2000 
April 20 1 1 :44:3 1^, and its gamma-ray light curve is shown in 
Fig.|2l Note that I did not detect, using a Scargle (1982) peri- 
odogram, any periodic signal in the light curve within ^ 100 s 
after the GRB spike. 

The second match, GRB 000526B coincides with the posi- 
tion of the galaxy NGC 7331. However, as seen in Figure [8] 
this galaxy has an elongated shape (while in the search I as- 
sumed the galaxies are round) and the IPN annulus does not 
coincide with the galaxy apparent disk or any star forming 
regions within this galaxy. I note there is a smaller galaxy 
within the NGC 7331 group that does fall within the IPN 
annulus. However, this galaxy does not pass the absolute 
magnitude cut, and the a-priori probability to find an SGR 
flare from such a galaxy is small. Therefore, I conclude that 
GRB 000526B is unUkely to be associated with NGC 7331. 

I have estimated the probability for a chance coincidence 
between the short-GRB IPN error regions and the galaxy sam- 
ple. This was done by randomizing the positions of the 316 
galaxies over the celestial sphere, and searching for a match 
with one of the IPN short-GRBs in the list. I repeated this 

* Note that error regions for GRB positions obtained by some of the space- 
craft are not circular 

' On April 20th, M74 is about six deg from the Sun 
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Fig. 6. — A GALEX near UV-band image of the Sc-type galaxy M74 
(NGC 628). The parallel lines mark the section of the IPN annulus of 
GRB 000420B that passes through the galaxy. The IPN constraints for 
this burst include an annulus whose parameters are: annulus center, a = 
303.76 deg, S = -29.4 deg (J2000.0), annulus radius of 89.63 deg, and semi- 
width of 0.04 deg; and a KonusAVind ecliptic latitude constraints which indi- 
cates that the GRB ecliptic latitude is between +2 deg to +22 deg (3-0'; The 
ecliptic latitude of M74 is +5.2 deg). The total area of this annulus segment 
is 3.25 deg-. 
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Fig. 7 . — The KonusA¥ind gamma-ray light curve of GRB 000420B, which 
might have originated in M74. The GRB shows a single, short (0.3 s width), 
peak. 

simulation 1000 times and found that the expectancy number 
of matches is 3 and the probability to find two matches is 
^ 15%, and three or more matches is ^ 76%. Note that since 
the actual galaxy positions are correlated and the IPN con- 
straints are coordinate dependent (e.g., Konus/Wind ecliptic 
latitude constraint) this simulation does not mimic the search 
process accurately. However, it indicates that such a chance 
coincidence is possible. Therefore, I cannot securely identify 
GRB 000420B with M74. 

Given the product of the completeness factors, EqEsfr^-v ~ 
64% (see ij3j, and conservatively assuming that the search 




Fig. 8.— R-band image of the galaxy NGC 7331 from the Palomar Ob- 
servatory Sky Survey II. The parallel solid lines mark the section of the IPN 
annulus of GRB 000526B that passes near this galaxy. The center of the 
dashed-line circle around NGC 733 1 shows the actual position of NGC 733 1 
in the Tully (1988) catalog, which is accurate only to about one arcminute. To 
avoid any incompleteness caused by the positional inaccuracies, I increased 
the galaxy radii by one arcminute compared to the radii listed in Tully (1988). 
The circle size corresponds to the increased-galaxy radius. 

yielded a single SGR candidate among 47 short GRBs, I put 
an upper limit on the fraction of SGRs among short GRBs of 
< 16% (< 27%), at the one sided 95% (99.73%) CL (Gehrels 
1986). Assuming that GRB 000420B is indeed an SGR flare 
in M74, the fraction of SGRs among short GRBs is about 3% 
\ ) 

5. THE RATE OF GIANT FLARES 

A closely related quantity to the fraction of SGR giant flares 
among short GRBs is the rate of SGR giant flares, which I 
estimate below. 

Among the four active SGRs in the Galaxy (including the 
LMC), a single giant flare with energy above 3.7 x 10^^ erg 
was observed in the last 30 years. Therefore, the derived 
rate of giant flares with energy above 3.7 x Vd^^r\c^ erg (Hur- 
ley et al. 2005) is (120t^^|j°)-' yr"' per SGR (95% confi- 
dence interval; CI; Gehrels 1986), where rjs is the distance to 
SGR 1806-20 in 15 kpc units (Corbel & EikenbeiTy 2004). 

Next, using the results of 2] I estimate the rate of SGR 
giant flares, and compare it with the rate derived from the 
Galactic SGR giant flare. Assuming that all the SGRs have 
identical flare properties (i.e., the same cut-off energy Zimaj 
and a), then the rate of giant flares, with energy above is^/, 
per year per SGR, /g/, is given by: 

(3) 

where Ar is the time span of observations (15 yrs), No,i is the 
observed number of SGR flares in the /-th galaxy (E.iNo.i < 1), 
a = 5/3 (Cheng et al. 1996), A^^'^'^ is the number of active 
SGRs in the /-th galaxy, and C is the completeness function 
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(Eq.[T]i. To estimate A^^*'^'^, I use the total far-IR flux (see fj3j 
of the galaxies (Table[lli as an estimator for the star formation 
rate (Helou et al. 1988). I converted the far-IR flux to the SN 
rate, Rf', by normalizing the far-IR flux of each galaxy to that 
of NGC 253, which I assumed has 0.2 SN yr"' (Pietsch et al. 
2001). The far-IR flux along with the derived SN rate for each 
galaxy in this sample are listed in Table[T] Finally, the number 
of SGRs in the i-th galaxy, A^"'"''^, is: 

A^r=A^«;:— , (4) 

where, Nmw is the number of active SGRs in the Milky-Way 
galaxy (four), and the SN rate in our Galaxy, R'™„,, was set to 
0.03 yr-i. 

Now I can compare the observed Galactic rate of giant 
flares with the rate implied by the extragalactic survey. Fig- 
ure |9] shows contours of equal fgf/Y.iNo,i (i.e., the rate fgf 
in Eq.|3]calculated assuming a single SGR flare was detected 
in the sample, £, A^„., = 1), as a function of the giant flare 
energy Egf and the cut-off energy, E,„ax- The labels on the 
contours denote the SGR giant flare rate in yr^' per SGR. 
The 95% confidence region on the rate of giant flares, with 
energy above 3.7 x 10^^ erg, as implied by the Galactic SGR 
observation (i.e., 0.0083^ooQg[ yr^^ per SGR) is shown as 
the light-gray region on the left side of the Figure, and the 2- 
(7 lower limit (i.e., approximately 2 x lO^'* yr^') is marked as 
a bold line. The estimated energy with its uncertainty, of the 
2004 December 27 giant flare, assuming a distance of 15 kpc 
to SGR 1806—20, is marked as dark-gray region on the right 
side of the figure. For convenience, I plot on the upper axis of 
Figure|9]the distance to SGR 1806—20 which corresponds to 
the energy of the 2004 December 27 giant flare. Plotted on the 
right-axis is the approximate magnetic field, B, corresponding 

to E„,ax i— §^^'^RIs^ where I set the neutron star radius, R^s, 
to ten km). 

Apparently, there is a conflict between the Galactic rate of 
SGR giant flares and the rate implied by the extragalactic sur- 
vey which predicts a giant flare rate of ^ 0.5 x 10^^ yr^' (as- 
suming Egf ~ 3.7 X 10*^ erg). There are several ways to ex- 
plain this contradiction: The observed number of SGR flares 
in the search, Y^jNoj = 1, is subject to Poisson error llggg 
(95% CL); The estimate of N'-^' (Eq. |4|i is subject to errors 
through the SN rate estimate, and to Poisson error in the num- 
ber of active SGRs in the Galaxy (4+f ; 95% CL); The dis- 
tance to SGR 1806—20 may be smaller than 15 kpc (e.g., 
12kpc, Figer et al. 2004; 6-10 kpc, Cameron et al. 2005), 
and therefore shifting the dark-gray zone (i.e., energy of the 
2004 December 27 giant flare) in Fig.|9]to the left; Moreover, 
the value of the power-law index a may be different from 5/3 
(e.g., Gotz et al. 2006; a ~ 1.9), and the luminosity function 
model is an approximation. In reality the energy distribution 
(i.e., dN/dE °^ E^^l^) is probably not terminated abruptly at 
E^ax, as different SGRs may be born with different magnetic 
fields. I note that although the tentative detection of a single 
extragalactic SGR has a large Poisson uncertainty (l^g gg, at 
the 95% CL), the actual rate cannot be much smaller than one, 
or it will not be consistent with the observed Galactic rate of 
giant flares with energy above 3.7 x 10^^ erg. 

Interestingly, the Galactic rate of giant flares with energy 

above 2 x lO'*'^ erg is (40+20-2^)"' yr"' Per SGR (The er- 



rors are the 1 and 2-a CL, respectively; assuming three^ gi- 
ant flares were observed in 30 yrs). This rate is indicated by 
the error-bar in the left-hand side of Fig. |9l The solid line is 
the l-(7 CI, and the dashed line is the 2-a CI implied by the 
Galactic rate. In order to reconcile the Galactic rate with the 
extragalactic rate, a gradual break (or steepening) in the flare 
energy distribution is required at < 3 x 10^* erg (at the 95% 
CL). 

I note that if I adopt the value of a = 1 .9 found by Gotz et 
al. (2006), then the inferred rate of giant flares with energy 
above 4 x lO'*^ erg will be even lower (by about 10%) than 
the rate based on a = 5/3. Moreover, for o: = 1.9 I find that 
Emax < 3 X 10"*^ at the 95% CL. 

To conclude, this analysis shows that: E,nax Si 3 x 10^^ erg 
at the 95% CL; the rate of SGR giant flares is (0.4 - 5) x 
lO^'* yr^' per SGR; and the fraction of SGR flares among 
short-duration GRBs cannot be much smaller than one per- 
cent, otherwise it will be inconsistent with the Galactic rate. 

6. DISCUSSION 

The fraction of SGR flares among short GRBs and the rate 
of SGR giant flares are closely related issues. Measurements 
of these quantities are important for the understanding of the 
SGR flare mechanism and magnetic field strength. Below I 
discuss the possibility of observing additional extragalactic 
flares and the implications of the measured SGR giant flare 
rate for the magnetic field strength in the context of the mag- 
netar model. 

The derived fraction of SGRs among short GRBs implies 
that about five to ~ 65 extragalactic SGR flares were ob- 
served by BATSE. Moreover, I estimate that the Swift satel- 
lite may find an extragalactic SGR every several years, and 
that the Gamma-Ray Large Area Space Telescope will lo- 
calize about 0.2-3 extragalactic SGRs per year Based on 
the analysis presented in this paper, it would be surprising if 
both GRB 000420B and the recently discovered GRB 05 1 103 
(Golenetskii et al. 2005; Ofek et al. 2006; Frederiks et al. 
2006) were not associated with M74 and M81, respectively. 

Verifying that a GRB is indeed an extragalactic SGR flare 
requires either the detection of persistent X-ray emission, 
which is beyond our current observational capabilities, or 
the detection of additional flare activity. The fact that SGR 
flares are correlated may increase the probability of observ- 
ing additional flare activity from extragalactic SGRs. The 
observed time between low-energy SGR flares follows a log- 
normal distribution (Laros et al. 1987; Hurley et al. 1994; 
Cheng et al. 1996). Therefore, it is more appropriate to ex- 
press the typical time between bursts using the log-normal- 
mean, M, and -standard deviation, S. Assuming a flare rate 
of 2 X 10^^ yr^^ then T, the expectation value of the time 
between giant flares, with energy above 3.7 x 10*^ erg, is 
5000 yr. Since the expectancy value of log-normal distri- 
bution is given by exp(M + 5^/2), and assuming S = 3.46 
(Hurley et al. 1994) is independent of energy, the log mean 
for giant flares is Mgj = 2.53. Interestingly, the cumulative 
log-normal distribution suggests that given these parameters 
(i.e., Mgf = 2.53; S = 3.46), the probability to observe a 
second giant flare of an SGR within ten years from the first 
flare (given a sufficient energy source) is 47%. Therefore, if 

^ These are: 1979 March 5 (SGR 0526-66); 1998 August 27 
(SGR 1900+14); and 2004 December 27 (SGR 1806-20). Note that, SGR 
flares are correlated (see also f|6) and therefore Poisson statistics can be used 
only if the flares were originated from different objects - as in this case. 
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Fig. 9. — Contours of equal /g//L, A'o.i, as a function of the giant flare energy, Egf, and the cut-off energy, Ei„a_x. The labels on the contours are the SGR 
giant flares rate in yr"' (i.e., rate of events for which E > Egf). The triangular area in the bottom-light is the unphysical region in which Egf > E„,ax- The 95% 
confidence region of giant flare rate (per SGR), with E > 3.7 x lO''* erg, as implied by the Galactic SGRs observations, (i.e., 0.0083^jJ jj^gj yr^'), is shown as 
the light-gray region on the left side of the Figure, and the 2-fT lower limit (i.e., approximately 2 x 10"'' yr"') is marked by a bold line. The estimated energy of 
the 2004 December 27 giant flare, (3.7 ±0.9) x lO"*^ erg assuming a distance of 15 kpc to SGR 1806—20, is marked as a dark-gray region on the right side of 
the Figure. The error-bar on the left-hand side marks the rate of Galactic SGR giant flares, with E >2x erg. The l-fj CI is marked by a solid line and the 
2-cT CI is marked by a dashed-line. For convenience, I plot on the upper axis of the Figure the distance to SGR 1806—20, which corresponds to the energy Egf, 

of the 2004 December 27 giant flare. On the right-axis I plot the approximate magnetic field, B, con'esponding to E,,,,,-, (= ^ 5^^ns)> where I set the neutron star 
radius, R„,, to ten km. 



GRB 000420B is associated with M74, and/or GRB 051103 
is associated with M81, there is a fair chance to see another 
flare from these galaxies in the next decade. Furthermore, 
I note that if more than one burst per SGR will be observed, 
the use of Poisson statistics will considerably overestimate the 
flare rate. 

As discussed in iJT] Stella et al. (2005) ai-gued that about 70 
giant flares, with energy above 5 x 10"*^ erg, are expected from 
an SGR during its life-time, and claimed that SGR magnetic 
fields are in excess of 10^^ G. My analysis, however, shows 
that the giant flare rate is two orders of magnitude smaller than 
that estimated by Stella et al. (2005), and therefore, magnetic 
fields of B « 10^^ G can provide the required energy source. 
The giant flare rate I derive, (0.4 - 5) x 10""^ yr"' per SGR, 
suggests that within the lifetime of an SGR, 10^-10^ yrs, it 
would have of order unity giant flares with energy above 4 x 
10"*^ erg. 

To summarize, I searched for a spatial coincidence of IPN 
GRBs with galaxies within 20 Mpc. I have found a single 
SGR candidate, possibly located within M74. However, this 
match could be a chance coincidence. Currently, the best way 
to test the hypothesis that GRB 000420B arose in M74 is to 
detect additional SGR flares from this galaxy. As SGR activ- 
ity is correlated, the expected waiting time to a second flare 
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